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1 Introduction 

Extensions of the SM often include one or more additional abelian gauge symmetries beside 
the Standard Model U(l)y hypercharge. Compelling motivation for these extensions arises 
from grand unified theories (GUT) and from string theory. The phenomenology of new 
abelian gauge groups has been widely studied in the literature [1-11]. 

Moreover, in the recent years there was a renewed interest in this kind of models, based 
on the key observation [12, 13] that an extra U(1)d GeV gauge boson would be able to 
naturally explain the anomalies observed in indirect dark matter detection experiments, like 
the reported PAMELA result on the positron fraction [14]. The extra gauge symmetry may 
be hidden from the SM particles, which are singlets under the new force: they constitute 
the visible sector. The particles charged under the new U(l)u and singlets under the SM 
gauge groups are often referred to as the dark (or hidden, or secluded) sector. 

The connection between the dark and visible sectors is established through mixing 
operators. One candidate term is kinetic mixing of U(1)y with U(1)d [15, 16]. Since 
cosmological considerations (like Big Bang Nucleosynthesis) severely constrain a massless 
gauge boson, the extra U(l) symmetry has to be broken. Its breakdown can be achieved 
through the introduction of a new Higgs boson, ho, which can naturally mix with the SM 
Higgs [17-21], thus providing an extra portal between the two worlds. These extensions 
of the SM have also been studied in the context of electroweak phase transition (EWPT) 
[22-24] and dark matter (DM) since the dark sector provides natural DM candidates [25— 



While the extra U{\) models are very constrained from current experimental data 
[46-50] , the non-observation of a Higgs boson yields very mild bounds on the Higgs portal 



45]. 
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parameters at present. In this paper we explore the constraints and detectability prospects 
of the Higgs sector at colliders. Very recent work [51] was focused in the potential signatures 
at the LHC for large luminosities (O(30) fb _1 , see also Refs. [52-59] for older studies). Our 
interest resides in the reach of the early LHC data (-^/s = 7 TeV, with a total integrated 
luminosity less or equal than 15 fb _1 ). Similar work was already done in the context of 
the MSSM in ref [60]. 

Due to the mixing with hypercharge, the dark gauge boson can be in conflict with 
electroweak precision data, such as the Z mass or the effective weak mixing angle. Thus 
the most natural options, already considered in the literature, is to have either a very 
heavy (TeV scale) Z' [61-63] or very light (GeV) boson [64-67]. The latter scenario is well- 
motivated when looking to find a unified explanation of recent results of DM, as suggested 
in ref [12]. In this work we take an agnostic attitude and consider the Z' mass a free 
parameter. 

This paper is organized as follows: in Section 2 we review the model under considera- 
tion. In Section 3 we explain in detail the scan of the parameter space, all the constraints 
under consideration and the LHC expected reach for different scenarios. Section 4 contains 
the numerical results of our analysis. Finally, we conclude in Section 5. 

2 Model review 

Generic dark sector models were discussed in detail in the literature (e.g. [19]). In this 
section we briefly review the model used in our study. The Lagrangian can be written as 
follows, 

£ = C-SM + C-Dark + fornix > 

(2.1) 

where we have split the contribution into the SM-piece, the dark sector and the mixing 
between the two sectors. For the dark sector, we would like to add the minimum field 
content. Thus, we include a new dark gauge boson X and a dark Higgs field Ho- The dark 
sector might contain fermions, which are SM singlets and charged under U(1)d- These 
fermions are, however, irrelevant in the present context. The dark Higgs field will give 
mass to the X boson after spontaneous breakdown of the gauge symmetry. We pick a U(l) 
gauge group for simplicity; that is not to say that the dark sector has to be that simple, 
but that we choose to parametrize it in a simple way. It is clear that many other, richer 
possibilities (from a phenomenological point of view) can also be considered 1 . Under these 
assumptions, the dark Lagrangian reads 

C D ark = (D^H D )\D^H D ) + ^ D H ] D H D - \ D {H D H ] D ) 2 - ^X^X^ + ... (2.2) 

where = + ig'YB^ + igT a W® + ignQoX^ is the covariant derivative, go the dark 
U(1)d gauge coupling, X^ v its gauge strength tensor and Qd is the charge under the dark 

1 One could argue that the details of the dark sector at energies above LEP and SLC could be absorbed 
into the low energy GeV scale parameters by integrating out the heavy sector. Another option is to work 
with a different dark gauge group. We will stick, for the sake of simplicity, to this minimum extra added 
field content. 
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force. The last term is the kinetic term for the dark gauge field, while the remaining terms 
correspond to the kinetic term for the complex scalar Higgs and the dark Higgs potential. 
The ellipsis stands for other terms not relevant for our study. The mixed Lagrangian will 
depend upon how X and Hp couple to the SM. In our setup, it is natural to consider 
kinetic mixing between X^ and and a mixing term in the Higgs potential [17, 18], since 
these two are the only renormalizable operators relating X and Hp to the SM. 2 With these 
assumptions we have 

Cmix = y B ^ X ^ + e H (HH^)(H D Hl), (2.3) 

where B^ v is the U(l)y hypercharge field strength tensor and H is the SM Higgs doublet. 
It is well known that has to be small in order to be compatible with current experimental 
limits (see [49] and references therein). The constraints on en are less stringent, given our 
current knowledge of the Higgs sector. 

For the sake of completeness, we write down the SM Lagrangian, 

C S m = {D^(D^H) + nHrf - A(Mt)2 

+ Y^Vf^HfR + h.c) - \{B^ + W£W«) + . . . , (2.4) 
/ 

where y/ is the Yukawa coupling for the SM fermion / and the ellipsis indicates the presence 
of other terms not relevant for our study. 

2.1 Gauge sector 

In order to derive the interactions in the mass eigenstate basis, we have to proceed in 
several steps. First, one has to diagonalize the kinetic terms for the gauge bosons. This 
can be achieved by performing a field redefinition of B^ and X^. After this, one finds 
that the covariant derivative has changed in such a way that now the dark sector interacts 
directly with the B^. Since we want the U(1)d gauge group to be broken, the vacuum 
expectation value of Bd will contribute to the masses of the Z and the Z', while the photon 
will remain massless. 

The Lagrangian involving both U(l) strength tensors is given by 

1 

4 1 

In order to diagonalize the kinetic term, we perform the following redefinition of the fields 
[11] first: 

Bfj, Bp + €A X M , X M -> 1 X^ . (2.6) 

Then the covariant derivative reads 3 

Dp = 8^ + ig'YBp + igT A Wp + i [ g D + g' | X M , (2.7) 



C = -- (B^B^ + X^X"" - 2e A B^Xp U ) . (2.5) 




2 As noted in ref [57], there are other such operators if the dark fermions are also taken into account. 
3 In our convention, the SM Higgs doublet has a Y—+1/2, and the dark Higgs doublet also has Qn = 
+ 1/2. 
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and the mass matrix of the neutral gauge bosons becomes 

/ 4 -cwsw asyy \ 
m z 1 -cwsw <?w -ac w s w I , (2.8) 

\ aSyy —CICWSW ^Syy + A' 

where sw , cw are the sine and cosine of the usual SM electroweak mixing angle and 

2 2 2 

2 / 2 , I2\ v 2 2 V D a mx e A / n \ 

m *> = (9+9) T , rn Xo = 9DW -^ y A = - | -, a = . (2.9) 



One of the mass eigenvalues is zero, corresponding to the photon eigenstate, and the two 
others are given by 

2 

M 2 = ^ [(1 + s 2 w a 2 + A) ± ^(l + ^a 2 + A)2-4A] . (2.10) 

Due to the smallness of ca it is well justified to take the gauge boson masses at their 
tree level values, namely, to assume mz = niz and mz> = mx - We have numerically 
checked that this approximation has an error below 0.02 %. The relation between mass 
and interaction eigenstates is given by 



/c w ~swc x s w s x \ (A 



] = ( s w c w c x -c w s x ] [ Z M ] , (2.11) 
and the new gauge boson mixing angle by 



' \ s x c x / \ Z'^ 



2.2 Higgs sector 

In the unitary gauge, one has 



B = Mvlh)' *> = >°+M. (2-13) 



and the minimization of the Higgs potential yields 



H = \v 2 -e H -^-, HD = ^Dvj) - e#y . (2.14) 



The squared mass matrix of the Higgs sector reads 



M 2 =( 2 ^ -fTV (2-15) 
y-e H vv D 2\ D v z D J 



with its eigenvalues given by 



m\ 2 = Xv 2 + \ D v 2 D T a/( Aw2 " A ^d) 2 + . (2-16) 
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where mi > mi. The mass eigenstates read 

h 2 = c a /i — s a h,D , hi = s a h + c a ho 

and the mixing angle is given by 
€ H vdv 



S2o 



(Xv 2 - \ D v 2 D ) 2 + e 2 H v 2 v 2 D 



C2c 



At; 2 - X D v 2 D 



(Xv 2 - X D v 2 D ) 2 + e 2 H v 2 v 2 D 



(2.17) 



(2.18) 



We define the effective Higgs coupling as the coupling in our model normalized to the SM 
case. Using eq. (2.17) in eqs. (2.4) and (2.2), one has 



9hiWW = 9h x j] = s » i 9h 2 WW = 9h 2 ff = c " 

The couplings to Z — Z' read 



(2.19) 



v v 
9h 2 z x z 2 = Caghz^ - Sa9h D z 1 z 2 A — ; 9h 1 z 1 z 2 = s a ghz 1 z 2 + c a gh D z 1 z 2 A — , (2.20) 

''/) " " vd 

where Z\^ = Z,Z' , the guz 1 z 2 factors are given in table 1. Due to the smallness of the 
kinetic mixing one finds that gnzz ~ gh D z'Z' ~ 1 ; while all the other are at least suppressed 
by a power of < 0.03. Therefore, one has that the coupling of h\ (/12) to the SM particles 
is suppressed by a factor of c a (s a ) with respect to the values of the SM Higgs. 



H 


ZZ 


Z'Z' 


ZZ' 


h 


(-c x + as x s w ) 2 


(s x + ac x s w ) 2 


(-c x + as x s w )(s x + ac x s w ) 


hd 


s 2 


c 2 


s x c x 



Table 1. gHZ ± z 2 couplings. 



There are also interactions involving three and four Higgs fields, as well as two gauge 
bosons plus two Higgs fields. These decay modes constitute what we will call, from now 
on, non-standard (Non-SM) Higgs decay modes, namely, those that do not appear when 
considering the SM Higgs boson. They could be important, for instance, if there is a 
significant fraction in the /12 — > h\h\ or /12 — > Z'Z' at LEP, like in the buried Higgs 
scenario [68]. In our setup we assumed that the decay width of the Z' into Standard Model 
particles is negligible, since its couplings to Standard Model particles are suppressed by a 
factor of €a with respect to the couplings of the Z. The decay width of a Higgs boson into 
two gauge bosons Z\ and Z2 is given by 



T(H -> Z X Z 2 ) 



9 m Hg H z 1 z 2 s 



rrir 



1 1 
m Zi m z 2 



[x\ + x 2 ) (x\ — Xl 



+ 



1 + -X\X2 + 



i6~~ 



4 

X\ + x 2 



1/2 



(2.21) 



where H = hi,h,2, £1,2 = (2?riZi 2 /m//) 2 , gHZ^z 2 can be read from table 1 and S is a 
symmetry factor, 1/2 if Z\ = Z2, 1 otherwise. The partial widths of the heavy Higgs boson 
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into light ones is 



T(h 2 



Am 2 hi 



32irmh 2 



in 



9fi2hihi 



(2.22) 



hi 



where the trilinear Higgs coupling gh 2 h 1 h 1 is given by 



2{3s a c a (Xvs a - XdVdCo) ~ x [vc a (3c 2a - 1) + v D s a (3c 2a + !)]}• (2.23) 



Due to the rescaling of the Higgs-to-Standard Model couplings the Higgs production cross 
sections are suppressed by a factor of s& for hi (c„ for h 2 ). Consequently, in the case 
where one can neglect the non-SM decays, there is always one Higgs boson whose rate is 
suppressed at most by a factor of 1/2. The branching fractions into SM particles will be 
suppressed by a factor of 1 — Br(/ij — > non-SM). Therefore, the total rate for any Higgs 
boson into SM particles is always lower than in the SM by a factor of 



sLww (1 ~ Br ( h i -> non-SM)) . 



(2.24) 



3 Numerical analysis: parameter scans and constraints 
3.1 Parameter scans and pre-LHC constraints 

To explore the parameter space of the model a random parameter scan was performed 
using the CliBA-library [69]. We chose as input parameters the physical parameters mi, 
m 2 , the mixing angle a, go, m,z' an d the kinetic mixing parameter €a with values in the 
ranges according to table 2. We focused on Higgs masses below 600 GeV since the LHC 



772! [GeV] 


m 2 [GeV] 


a 


m z > [GeV] 


9D 




[1;400] 


[1;600] 


[0;tt] 


[0; 1000] 


[0;i] 


[0;0.3] 



Table 2. Ranges of the parameter scan. 

experiments have published exclusions in that mass range and the phenomenology of a 
heavy singlet Higgs has been studied elsewhere (see, for instance, ref [54]). 
The potential parameters were computed using 



A 



1 

4v 2 

1 



[mj(l - c 2a ) +m|(l + c 2a )] , 
[ml (1 + c 2a ) + m|(l - c 2a )] 



2vvd 



{ml - m\) s 2a . 



(3.1) 
(3.2) 
(3.3) 



We also required the points to respect the positivity conditions, eq. (2.14), thereby ensur- 
ing the proper minimalization of the potential. Motivated by the discussion of the elec- 
troweak phase transition in similar models (see, for example, [70] and references therein), 
we discarded points with nonperturbative potential parameters by requiring en ^ 0.5 and 



- 6 - 



A, Ad ^ 1. This also limits the contribution of the invisible decay modes to the total width 
of the Higgs bosons such that their values stay within the validity of the narrow width 
approximation (i.e. r* ot /m.j < .05), which is required in order to interpret the exclusion 
limits set by collider data on the rates of the Higgs boson as the product of the production 
cross-section times branching ratio in a particular channel. 

Constraints from direct searches were applied using HiggsBounds 2 . 1 . 1 [71, 72], where 
points are excluded at the 95% confidence level. In the low-mass region (below 114.4 GeV) 
the main exclusion channels are the LEP searches for a Standard Model-like Higgs [73-75] 
and a Higgs- like scalar decaying completely invisibly [76-78]. In some cases the decay 
/i2 —> h\h\ —> 46 or 4r was also constrained directly by the corresponding LEP MSSM 
searches [73]. In the high-mass region (120-200 GeV) the Tevatron searches were also used 
to bound the parameter space [79, 80]. 

Electroweak precision data also limit the parameter space of our model in a significant 
way [81, 82]. In order to assess the effect of a complete parameter fit, we used model 
independent bounds on the kinetic Z — Z' - mixing [49] to constrain e A and computed 
the contribution of the extended gauge and Higgs sectors to the Peskin-Takeuchi S and T 
parameters [83] using FormCalc [84]. For the two Higgs bosons hi and h2, it is given by 

S = clS SM (mi) + slS SM (m 2 ), (3.4) 

where S SM denotes the contribution of a Standard Model Higgs with respect to the reference 
mass m/j = 120 GeV (and analogously for the T-parameter). The tree- level contribution 
of the Z' to the oblique parameters is [89] 

olewS = 4cwSw -c^ e A , and a E wT = -sfy A e A , (3.5) 

which diverge as mz* — > mz- We are however confident that the formulae are valid as 
long as It^iaH <C 1. Due to the constraints on e A that we implemented this condition is 
always fulfilled in our scan. Since we study the Higgs sector of this theory at the LHC, 
we were interested in how a Z' with suitably chosen properties can relax the upper mass 
limit on the Standard Model Higgs mass from the S and T parameter fit, which is the 
case when the tree level contributions are enhanced through mz 1 — > fnz- We neglect loop 
contributions of the Z' via the ordinary photon and W,Z gauge boson self energies, since 
their size would be of order of the Standard Model gauge sector contributions to the neutral 
current amplitude, but suppressed by an additional factor of e A . In the threshold region 
around mz> = mz it is suppressed even further by the strong constraints on e A . We also 
neglected dark fermions, since their contribution would only enter the forementioned at 
the two-loop level. We set U = and required a parameter space point to lie inside the 2a 
contour in the S — T-plane provided by the Gfitter collaboration [85-88] . 

3.2 LHC data and projections 

In our analysis we include the current LHC data and future projections for the search 
channels listed in table 3. All of the searches use the most up-to-date LHC data with a 
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Lum. (ftr 1 ) 
ATLAS CMS 


a t wp n o 

VV lldlj WO V-LVJ 


Mass range 
(GeV) 


Ref. 


pp^H -»• WW 


1.7 1.5 


Comb. 


115-600 


[90-92] 


pp-> H ZZ 


1.04-2.28 1.1-1.7 


Comb. 


120-600 


[93-99] 


pp — > H — > 77 


1.08 1.7 


Comb. 


110-150 


[100, 101] 


— >■ -ff — )■ t + t~ 


1.06 1.6 


Comb. 


100-150 


[102, 103] 


ViZ", H ^ bb 


1.1 


CMS x 2 


110-135 


[104] 


qqH, H — > t + t~ 


1 


ATLAS x 2 


110-130 


[105] 



Table 3. List of LHC channels used in this study. Here, H stands for either hi or /12. The production 
mechanisms considered in pp are gluon-fusion, vector boson fusion, associated production with 
Z, W,tt and also bb — > H. The cross sections at the LHC have been taken from ref [106]. Sec the 
main text for details. 

total integrated luminosity between 1.04 - 2.28 fb -1 (depending on the search channel), 
except for the qqH, H — > t + t~ channel, for which we use the MonteCarlo 2010 sample 
[105], which provides a projection of the expected sensitivity of the current data-sample, 
since no LHC collaboration has presented yet updated data in this search channel. In the 
case of the associated production with a vector boson, with the Higgs decaying into bottom 
pairs, the current analysis was done using a cut based procedure that is able to exclude a 
Higgs boson with a rate of around 20 times the SM case [107]. The MC 2010 analysis was 
performed by taking advantage of boosted bb pairs [108], and the expected exclusion with 
1 fb -1 of data for this case is around 6 times the SM [105]. We note that these channels 
are not able to probe points in our model, but, for the sake of completeness, we include 
them in our analysis. 

We combine the results from ATLAS and CMS in the channels where both collab- 
orations have presented data, following the prescription detailed in Refs. [109, 110] (see 
below). The current exclusion is obtained by using the observed limits reported by AT- 
LAS and CMS. For channels where only one of the collaborations has presented data, we 
will base our current exclusion on that analysis. In the same case, we compute the future 
projections by doubling the expected result in an attempt to mimic the combination of 
both experiments and scaling the result by the expected total integrated luminosity using 
the prescription detailed below. Since the reach for SM-like Higgs bosons at CMS and 
ATLAS is similar, this approximation is expected to be reasonably accurate. Except in 
the ZZ channel, all of the others involve, for a particular mass range, one definite final 
state. For the ZZ we also combine in quadrature the results for the a) four leptons, b) 
two leptons two quarks and c) two leptons plus two neutrinos, and d) two leptons plus two 
taus (CMS only) final state searches 

Here we briefly review the procedure used to combine the experimental information. 
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LHC @ 7 TeV, 15 fb _1 (ATLAS+CMS) 



LHC @ 7 TeV, 15 lb -1 (ATLAS+CMS) 




M H [GeV] M H [GeV] 

(a) (b) 

Figure 1. LHC reach (a) and statistical significances (b) for the SM Higgs boson H with 15 fb , 
combining both experiments. The color coding is as follows: WW (blue), ZZ (orange), 77 (red), 
tt (black) , VH, H — > bb (green) and qqH, H —> t + t~ (purple). 



For each relevant channel we compute the following quantity: 

Q(A>) = j^rh ' ^ 

where R mo d is the rate in this particular channel coming from our model, R e *p is the 
exclusion limit at the 95% C.L., at a reference total integrated luminosity £rj- Eq. (3.6) is 
exactly the same definition used by HiggsBounds in order to set the 2a exclusions: if Q > 1 
the point is excluded at the 95% C.L.. We compute R exp combining the results for each 
LHC experiment in inverse quadrature (see Refs. [109, 110]). In ref [60] this procedure was 
found to be more conservative than the combination performed by the ATLAS collaboration 
by 10-20 % 4 . While the quantity R mod is a number that does not change, i? eX p scales with 
the luminosity as C~ 1 / 2 . Thus, defining R exp (£o) = Rq and Q$ = Rmod/Ro, one has that 




In order to derive these equations, one is neglecting all systematic effects and also assumes 
that in each particular channel B 3> S S> 1 holds, where B (S) are the number of back- 
ground (signal) events for a particular channel. With these simplifications the expected 
statistical significance a turns out to be a w 2 Q. 

As an illustration, we present in figure 1 the expected reach at the LHC and the 
statistical significance for the SM Higgs as a function of the Higgs mass in the channels 



4 While the first version of this manuscript was under consideration, ATLAS and CMS presented the 
combination of their datasets in [111]. We have compared their results against our naive combinations, 
finding that the expected values differ by at most 10%, while for the observed values the discrepancy ranges 
from 30 to 50 %, but in those cases our naive combination turns out to be a conservative. A similar 
comparison is also shown in figure 2 of ref [112], where the SM combination is confronted against the 
experimental result, also finding a similar accuracy. 
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Figure 2. Scanned points in the mj — m-i plane. In the left panel (a) we show a representative 
sample of excluded along with allowed points, while in the right panel (b) we show only allowed 
points. In the left panel, we show the points excluded by LEP (blue), Tevatron (green), S- and T- 
parameters (red, filled), ca (red, empty) and perturbativity constraints (orange). The black points 
are allowed by all current collider data. In the right panel the color coding varies according to the 
value of c 2 a . red (0 < c 2 a < 0.25), magenta (0.25 < c\ < 0.5), blue (0.5 < c 2 a < 0.75) and black 
(0.75 < 4 < 1). 



described in table 3, assuming a total integrated luminosity of 15 fb , which corresponds to 
the total integrated luminosity that can be collected by the end of 2012 if the instantaneous 
luminosity is kept at the current rate. 

As one can see, the WW channel is setting the most stringent exclusion in the 120 — 200 
GeV range. For larger masses, the H — > ZZ channel takes the leading role. For masses 
below 120 GeV one enters into the problematic range, where the WW channel becomes 
ineffective, and the diphoton requires O(10 fb -1 ) integrated luminosity in order to probe 
the SM Higgs. Moreover, as mentioned in Section 2.2, the suppression factor for one of 
the Higgs bosons is at most 1/2, unless there is a significant non-standard branching ratio. 
This means that in the case of maximal mixing without significant extra-SM decay modes, 
at least one Higgs boson can be tested at the 2 (5) a level if its mass lies in the 125-550 
(140-190) GeV range. 

4 Numerical analysis: results 

In this section we present the results of the parameter scan as defined in Section 3. In 
the left panel of figure 2 we study the impact of each experimental or theoretical bound 
on the parameter space of the model in the mi — mi plane. We plot points that are 
excluded by LEP (blue), Tevatron (green), S and T parameter (red, filled), ea (red, empty) 
and the requirement of perturbativity of the potential parameters (orange). It is clearly 
visible how the direct searches of LEP (blue) and Tevatron (green) constrain the region 
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mi < 114.4 GeV and m% ~ 160 — 170 GeV. The indirect bounds via the S and T parameters 
(red, filled) and the constraints on kinetic mixing of the neutral massive gauge bosons (red, 
empty) mostly affect regions where one or both Higgs bosons are heavier than 155 GeV. 
As can be seen from equation (3.1), the perturbativity requirement A < 1 places an upper 
bound mi < V%v ~ 350 GeV on the mass of the lighter Higgs boson (orange). However, 
points in that region tend to be excluded for other reasons before so that its most important 
effect is to prevent the decay width of either Higgs boson into Z'Z' from becoming large 
enough to invalidate the narrow width approximation. 

The black points evade all of the above constraints and are the focus of the study at 
hand. This subset is shown in the right panel of figure 2, where we have colored the points 
according to the rescaled squared coupling of hi to Standard Model particles, namely, 
to the particular value of c 2 a : red (0 < (? a < 0.25), magenta (0.25 < <? a < 0.5), blue 
(0.5 < 4 < 0.75) and black (0.75 <c 2 a <l). 

We focus first on the region where the mass of the ligher Higgs state hi is below 
the LEP limit of 114.4 GeV. Here, the hi coupling to Standard Model particles has to be 
significantly suppressed to avoid direct detection 5 . The heavier state hi can be lighter than 
114.4 GeV at the same time if hi — > hihi — )■ SM is by far the dominant decay channel, in 
which case it can evade the constraints coming from the LEP searches for invisible decays 
of a Higgs-like scalar 6 . The decay hi — > ZZ' is subdominant (branching fraction below 1%) 
in this region as well as in the whole parameter space. 

When mi < 114 GeV < mi < 155 GeV, the heavy Higgs state hi behaves largely 
as the Standard Model Higgs, except for possible non-standard decays. In turn, if mi is 
above 155 GeV, fine-tuning between the S and T parameter contributions of the extended 
bosonic sectors is needed in order to evade the constraints from electroweak precision tests. 
It is necessary that niz' < mz, and for a given point in the mi-m,2-plane the allowed range 
for mz> becomes smaller, the larger mi is. 

When mi is above 114 GeV, there are no direct constraints on the mixing angle. The 
indirect constraints via the S and T parameter force the heavy Higgs to mostly decouple 
from the Standard Model for masses mi > 150 GeV. Where the decay hi — > hihi is 
possible, the corresponding branching fraction is always smaller than 0.5. BR(hi — >• Z'Z') 
can take any value, whereas BR(hi — > Z'Z') decreases with growing mi. The region where 
both Higgs masses are larger than 155 GeV is, again, the result of mz< < mz- 

Let us now study how the LHC experiments will probe the parameter space with their 
main search channels for the Standard Model Higgs boson. Even though production rates 
and decay widths are never enhanced in this model, the parameter space is already probed 
efficiently. As it can be seen on the left panel of figure 3, the current dataset (1.04-2.28 

5 The parameter space points where mi < 12 GeV should be taken with care, since the LEP search in 
the hiZ, hi — > bb is cut-off at this value [73] and there are other low energy experiments that can probe 
this mass range more efficiently than the searches included in this study (see also [113] for an analysis of 
the LHC reach.) . The detailed analysis of this region is outside the scope of the present work. 

6 Recent studies based on jet-substructure techniques show that the /i2 — > hihi — > ij final states can be 
tested at the 5a level at the 14 TeV LHC with O(10 - 100 fir 1 of data [114-116], depending on the model 
under consideration, and on mi and m2. Due to these reasons we decide not to include those analysis in 
the present work. 
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Figure 3. (a) LHC 2er exclusions for different scenarios and (b) most sensitive channels in the 
TOi — plane. In the left panel, we show the points that have been excluded by the LHC with the 
current dataset (green), and also those that can be excluded after collecting 5 fb _1 (red) and 15 fb -1 
(blue) of data. The black points will still be allowed by LHC data. In the right panel we remove 
the points excluded by the LHC today. Color coding varies according to the most sensitive channel 
to that particular point and filled (empty) shapes correspond to h% (hi). The pp —> hi —> WW 
decay mode is shown in red, the diphoton channel in blue and the ZZ channel in green. The black 
points are not sensitive to the LHC search channels under consideration. 



fb 1 ), shown in green, is able to exclude a vast majority of points in the region of mi > 
140 GeV, mostly due to the hi -> WW/ZZ channels. If mi < 114 GeV, h 2 -> WW/ZZ 
is probing values of m<i above 135 GeV range. With 5 fb -1 (red points) one can exclude 
almost the complete region mi > 130 GeV except in the case that h\ has small couplings 
to SM particles and at the same time a large branching fraction into Z 1 Z 1 or other invisible 
particles. With 15 fb -1 of data (expected in 2012) the diphoton channel will start to probe 
points in the 110-130 GeV range (blue points in the left plot). The region where both 
h% and hi are below the LEP limit (black points), and where the main decay mode for 
hi is into h\h\, can not be tested with the channels used in this study, since the LHC 
collaborations have not presented dedicated searches for this kind of decays (one would 
typically look into bbbb, T + r~bb or even t + t~t + t~). 

In figure 4 we show the 5c discovery potential of the model in the mi — mi plane 
(left panel) and the rate suppression factor (2.24) for the most sensitive search channel as 
a function of the Higgs mass which is more sensitive for exclusion/discovery at the LHC 
(right panel). 

From the left panel we see that a discovery with 15 fb -1 is only possible if either 
mass is larger than about 130 GeV. We have explicitly checked that the case where LHC 
discovers both Higgs states in the early run is only possible if the masses are in the range 
130 < mi < 170 GeV and 130 < mi < 260 GeV and the mixing between the two states is 
sizable. 
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(a) (b) 

Figure 4. Left panel: 5a discovery potential with 5 fb^ 1 (green), 15 fb^ 1 (red), exclusion at 
15 fb _1 (blue) and parameter space points outside of early LHC reach (black). Filled (empty) shapes 
correspond to the discovered or excluded particle being h 2 (hi). Right panel: Rate suppression 
factor (2.24) at the LHC as a function of the mass of the Higgs more likely to be detected first. 



When a Higgs scalar with a lowered rate a ■ BR than the Standard Model expectation 
is detected at the LHC it is a priori impossible to decide which of the two mass eigenstates 
has been discovered using above searches. Furthermore, more involved studies are needed 
to determine whether the origin of the rate suppression is the mixing between the states 
or decays invisible to the specific search channel, e.g. hi — > Z'Z' or hi — > h\h\. From the 
right panel we see that after collecting 15 fb -1 of data, one can exclude a rate which is 0.6 
(0.05) times the SM rate for m h ~ 130 (160) GeV; for m h > 200 GeV this value is 0.2. If a 
Higgs in the mass range of 180 < moiscovery ^ 300 GeV is discovered with a moderate rate 
suppression factor a • BR/SM « 0.7, it is likely that the detected particle is the lighter 
mass eigenstate h\. Furthermore, the early discovery of a Higgs state with a mass larger 
than 155 GeV points toward m' z < mz- This is because it is likely that the first Higgs to 
be discovered is the one that couples more strongly to the Standard Model. If it is heavier 
than 155 GeV, a Z' gauge boson with specific properties is needed to reconcile such a high 
Higgs mass with electroweak precision data. 

5 Conclusions 

In this work we have studied in detail the constraints on dark Higgs models at the 7 
TeV LHC. In the scenario under consideration the usual SM Higgs boson (i.e the one 
responsible for electroweak symmetry breaking) mixes with a complex singlet that breaks 
an extra U(1)d gauge symmetry, which in turn mixes with the SM hypercharge U(1)y- 
The free parameters of this model are the masses of the two Higgs bosons mi and m,2, the 
cosine of the mixing angle c a between them, the mass of the additional gauge boson mz> 
with the kinetic mixing parameter ca and the coupling strength g^,. A parameter scan was 
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performed and the effect of theoretical and experimental bounds from direct searches at 
LEP and Tevatron and electroweak precision data was studied. 

The current LHC published analyses (with luminosities between 1.04 and 2.28 fb _1 ) 
are able to exclude masses above 140 GeV, unless there is a significant mixing between 
the two mass eigenstates, while with 5 fb _1 most of the points with the lightest Higgs 
mass above 130 GeV are excluded. Furthermore, we have found that the 7 TeV LHC with 
15 fb _1 of total integrated luminosity will be able to probe most of the parameter space 
at the 2<t level. The points that evade these constraints correspond to two cases. The first 
of them is when either Higgs boson lies in the 115 — 130 GeV range and there is some 
non-negligible mixing between the two mass eigenstates. For such masses the WW final 
state is not very useful, and if there is some mixing between h± and \i2 one can loosen the 
exclusion power of the 77 channel, which will rule out that mass range for the SM Higgs 
boson. However, we would like to point out that, according to projection on the SM Higgs 
that include combinations of the different low-mass sensitive channels, this region can in 
principle be accessed during the 7 TeV run. 

The second case takes place when the Higgs rates are diminished due to a sizable 
mixing between the mass eigenstastes, or if /12 —¥ h\h\ is kinematically open. For this 
region of parameter space one should perform a dedicated search of h% — > h\h\ by looking 
at final states like bbbb, T + r~bb or even t + t~t + t~ . 

We would like to remark that the effectiveness of these bounds gets looser if there is 
an important partial width of any Higgs boson into other (not specified in this work) dark 
sector degrees of freedom (like the dark matter candidate). As a consequence, our results 
can either be interpreted as valid in a completion of the model where the aforementioned 
channel is not relevant, or also as the largest exclusion coverage one can get in parameter 
space. 

We have also analyzed the possibility of discovering one or two Higgs bosons. We have 
found that, with 15 fb _1 one can discover one of the Higgs bosons of this scenario if their 
masses are larger than 130 GeV. If only one Higgs boson with a rate smaller than the SM 
Higgs is discovered, it would be impossible to tell a priori whether it is h\ or hi- Discovering 
both Higgs bosons will only happen if 777-2 £ [130, 260] GeV and mi G [130, 170] GeV. Such 
an observation would rule out some other models, like for instance the MSSM, since the 
lightest neutral Higgs boson can not have a mass well above 135 GeV. Finally, if the LHC 
should not see any Higgs signature after collecting 15 fb _1 of data then one can use that 
result to constrain the mixing between the two Higgs bosons and/or the invisible width. 
For the former case, a recent study using the H — > ZZ — > 41 lineshape was presented in 
ref [117]. 
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Note added: while this paper was in press, the LHC collaborations have shown updated 
results using a larger dataset (about 5 fb -1 of data). Their results hint on a signal of a 
Higgs boson with a mass of around 125 GeV with rates compatible with the SM. These 
results leave our conclusions unchanged. 
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